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Abstract 

Due to the limited statistics so far accumulated in the Higgs boson search at the LHC, the 
Higgs boson property has not yet been tightly constrained and it is still allowed for the Higgs 
boson to decay invisibly to dark matter with a sizable branching ratio. In this work, we perform a 
comparative study for the Higgs decay to neutralino dark matter by considering three different low 
energy SUSY models: the minimal supersymmetric standard model (MSSM), the next-to-minimal 
supersymmetric standard models (NMSSM) and the nearly minimal supersymmetric standard 
model (nMSSM). Under current experimental constraints at 2c level (including the muon g — 2 
and the dark matter relic density), we scan over the parameter space of each model. Then in the 
allowed parameter space we calculate the branching ratio of the SM-like Higgs decay to neutralino 
dark matter and examine its observability at the LHC by considering three production channels: 
the weak boson fusion VV — h, the associated production with a Z-boson pp > hZ + X or a pair 
of top quarks pp > htt + X. We find that in the MSSM such a decay is far below the detectable 
level; while in both the NMSSM and nMSSM the decay branching ratio can be large enough to be 
observable at the LHC. We conclude that at the LHC the interplay of detecting such an invisible 


decay and the visible di-photon decay may allow for a discrimination of different SUSY models. 


PACS numbers: 14.80.Da,14.80.Ly,12.60.Jv 


I. INTRODUCTION 


As a cornerstone of the standard model (SM) and also the last undiscovered piece, the 
Higgs boson has been intensively searched in collider experiments. The foregone colliders 
LEP II and Tevatron yielded null search results, setting a lower bound of 114.4 GeV on 
the Higgs mass |l| and excluding a Higgs boson with a mass around 2My B. respectively. 
The ongoing Large Hadron Collider (LHC) took over the Higgs-hunting task and recently 
reported its search results. Based on an integrated luminosity of 4.9 fb ! collected at 
vs = 7TeV, the two experimental groups at the LHC independently further narrowed down 
the Higgs mass region (at 9596 C.L. the CMS collaboration excluded 127-600 GeV while the 
ATLAS collaboration excluded 112.9-115.5 GeV, 131-238 GeV and 251-466 GeV) and both 
hinted to a Higgs boson around 125 GeV |3]. Such a finding has stimulated some theoretical 
studies for a Higgs boson near 125 GeV in low energy supersymmetry |4| and other models 
5]. 

Of course, if the LHC hint of a 125 GeV Higgs from the di-photon channel is confirmed 
in the future, it would severely constrain or exclude those new physics models in which some 
new exotic decay modes (such as decaying invisibly into dark matter) are open and the di- 
photon rate is suppressed. But so far the statistics at the LHC is too small to confirm such a 
Higgs, let alone the precision measurement of the Higgs decay branching ratios. Therefore, 
experimentally it is still allowed for the Higgs boson to decay exotically, such as invisibly to 
dark matter, with a sizable branching ratio +. 

Theoretically, the Higgs decay to dark matter can indeed occur in some new physics 
models, such as the gauge singlet extensions of the SM [7], the SM with a heavy fourth 
generation |8], the large extra dimension model |9], the technicolor model " the spon- 
taneously broken R-parity models and the non-linearly realized supersymmetric model 
and the MSSM with a singlet . In this work, we perform a comparative study for the 
Higgs decay to neutralino dark matter in low energy SUSY by considering three different 
models: the minimal supersymmetric standard model (MSSM) [14H16], the next-to-minimal 


supersymmetric standard models (NMSSM) "n and the nearly minimal supersymmetric 


1 In " the authors used the limited LHC statistics of the di-photon signal rates to set constraints on the 
invisible Higgs decay in the MSSM and found that the invisible Higgs decay branching ratio around 1096 
is allowed. 


standard model (nMSSM) Lid by). As will be shown, in both the NMSSM and nMSSM, the 
SM-like Higgs boson can decay to neutralino dark matter with a sizable branching ratio. 

In case that the Higgs boson decays to dark matter with a sizable branching ratio, de- 
tecting such a decay at the LHC will be important because in this case the conventional 
visible decays into yy, bb, r7, WW? and ZZ“) are often suppressed. Obviously, the main 
production channel via gluon fusion gg — h is not usable because it just gives missing en- 
ergy. It was found through Monte Carlo simulations that the production via vector boson 
fusion (VBF) pp — hqq' and the associated productions pp — hZ and pp — htt can offer 
the opportunity to detect the Higgs decay to dark matter bdbd So in this work we choose 
these three production channels to display the observability of Higgs decay to dark matter 
in low energy SUSY. 

Note that although in the literature the Higgs decay to nutralino dark matter has been 
discussed in some specific model like MSSM, it is necessary to give a revisit in low energy 
SUSY for two reasons: (i) Different SUSY models usually give rather different phenomenol- 
ogy and it is interesting to perform a comparative study for different models; (ii) We want 
to know in the SUSY parameter space allowed by current experiments whether or not the 
Higgs decay to nutralino dark matter is detectable at the LHC. If in some model this decay 
is found to be accessible at the LHC, we further want to know how large the parameter 
space can be covered by searching for such an invisible decay at the LHC. 

This work is organized as follows. In Sec. II, we briefly describe the three supersymmetric 
models. In Sec. III, through a scan over the parameter space we present the branching ratio 
of Higgs decay to neutralino dark matter and show the observability at the LHC. Finally, 


the conclusion is given in Sec. IV. 


I. THE SUSY MODELS 


In a renormalizable supersymmetric field theory, the interactions and masses of all par- 
ticles are determined by their gauge transformation properties and superpotential. The 
superpotential is a holomorphic function of chiral superfields ®; 5 (ó;, v; F;), with aa 


and F; being respectively the bosonic, fermionic and auxiliary fields, and takes a form |15 


E T5. oem l cmm 
W = Lid, + 5M" 9,0; + gy 99,0, (1) 


where the parameters L; should be of dimension [mass]? and is only allowed if Ê; is a gauge 
singlet. The mass matrix M" can only be non-zero when the supermultiplets Ê; and ©; are 
conjugates of each other under gauge transformation. And the massless coefficients y'/^ can 
only be non-zero when $,0,0, formed a gauge singlet. 

The MSSM is the most economized realization of supersymmetry in particle physics, 


which has two Higgs doublets H,, and Hy and its superpotential is given by |15 
WMSSM — W, + uÊ- Hu, (2) 
with Wp given by 
Wr = uY,Q- H, — dY4Q - Ha — eY.L - Ha. (3) 


The MSSM has the so-called u-problem, which can be solved in some extensions by intro- 
ducing a Higgs singlet $. Among these extensions the most popular ones are the NMSSM 
and nMSSM, whose superpotentials are ld 


WNMSSM — Wp + AS, - Ha ;$* (4) 
W^MSSM — Wp + ASH, ] Hi, + &pM28. (5) 


The scalar potential in the Lagrangian contains the so-called F-term, D-term and soft-term 


2 
VEUSY = Vp + Vp + Voofts (6) 


where 


OW 


Vr = F“F,, F* = -WŻ = ——, 
0$; 


G? g2 
Vp = —-(HaP - |Hul?)? SE QBGPLBSP — |Ha + Hal’); (8 
Vote = Mi, Hul? + mh, Hal + (BH, - Ha + h-c); (9 


K 
Van am |Ha? + Th, | Ha + sls]? + (AAS Hu - Ha + gAs + h.c.) (10 


) 
) 
: ) 
VaMSSM — ge LH + mal Hal? + MISI + (XAXSH,- Hat EsM2Ê h.e) (11) 


Here the parameter G is defined as G? = g? + g2 with gı and g> denoting respectively the 
U(1)y and SU(2); couplings. With the Higgs fields H,, Hy and S developing respectively 


the VEV v, vg and vs, they can be rewritten as 


m= Ht P ud g va dbs $ ips n" 
u eines > ffd — H- , = V2 ( ) 
V2 d 


In both the NMSSM and nMSSM we have five complex scalar fields or ten real scalar degrees 


of freedom, whose mass eigenstates are obtained as 


h uU a U 
" 8 : ee ] e des (13) 
= Si | a | = Pij ; = Ci 
2 Qa à j | Pa Ae j H+ 
ha Qs G Ps 


Here the three Goldstone bosons G? and G^ will be eaten by the weak gauge bosons Z and 


W= respectively. Then we have seven Higgs bosons, among which ^4, ho and hg are CP-even 


(with the convention mp, < mj, < mj,), a; and a» are CP-odd (ordered as ma, < Mas), 


and H- are the charged ones. 
In the NMSSM and nMSSM there are five neutralinos (x2), which are the mixture of bino 
(B), wino (W°), higgsino (H,,, Ha) and singlino (S): 


xi B 
x3 we 
xg | =N] Ay (14) 
xy Ha 
x3 5 


We assume the lightest neutralino is the lightest supersymmetric particle (LSP) and make 
up of the cosmic dark matter. 
For the purpose of our numerical analysis, we present the interactions of the Higgs bosons 


h; with bb, rr, WW and x9X9. In the the NMSSM, they are given by [28 


hab dS : asi (15) 
hati Te : angel (16) 
hWiw, : a (aS + hasiz) (17) 
had: AGa + Sal Sal) ~ VZS Mi 

-A (Salli — Silli]) + 2 (Si Mii — Sif) (18) 


where II = Nu Ni; + Ni; Ny. The corresponding couplings in the nMSSM can be obtained 
by setting & equal to zero (for MSSM, setting «K, A, $13 and N15 to zero). 
For the NMSSM, in the basis X? = (B, W^, H,, Ha, $), the tree-level neutralino mass 


matrix takes the form |171, |28 


v2 v2 
0 M» = ETÀ 0 
Mo = an EE 0 =u —Auq (19) 
NI a u 0  —Au, 
0 0 | —Avug —Avy 2KS 


The corresponding mass matrix for the MSSM can be obtained by taking the upper 4 x 4 
matrix from the above expression, and for the nMSSM it can be obtained by set the term 
2&S to zero. When |p.rr| or Mo > Mz, the lightest neutralino in the MSSM becomes 
bino-like, with a mass given by » 


m? sin? 6,,(M, + usin 28) 
p? — Mj 


In the nMSSM the lightest neutralino is singlino-like and its mass can be approximated as 


(20) 


2ur? (v2 + v2) tan 8 


xv —— 21 
Tg 2u? + M (v2 + v2) tan? B +1 (21) 


III. NUMERICAL RESULTS AND DISCUSSIONS 


We scan over the parameter space of each model under current experiment constraints, 
and for each survived sample we calculate the Higgs spectrum, decay branching ratios and 
production rates at the LHC. In our calculation we use the package NMSSM Tools and 
extend it to the nMSSM bi. For the calculation of h — ^y in the SM, we use the package 

" (this 
code is aimed at the MSSM, and we extend it to the NMSSM and the nMSSM). For parton 
distributions we use CTEQ6L with the renormalization scale and the factorization scale 


Hdecay |29]. For the Higgs production cross sections, we use the code on the website 


chosen to be the sum of the masses of the produced particles. 
In our scan we require the models to explain the cosmic dark matter relic density measured 
by WMAP |32] and also explain the muon anomalous magnetic moment at 2c level. In 


addition, we consider the following experimental constraints: 
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1. The LEP bounds on sparticle masses and on the Higgs sector from ete~ — hZ(hA) 
followed by Z > £*£7,x0x9, h — bb and r*7- 33. We also consider the LEP- 
I constraints on the invisible Z decay ?, i.e., T(Z — x9x9) < 1.76 MeV, and the 
LEP-II constraints on neutralino production a(ete~ > x9x?) < 107? pb (i > 1) and 


a(ete” — x9x9?) < 107! pb (i,j » 1 33. 
j 


i 


2. The Tevatron bounds on sparticle masses and on stop or sbottom pair production 


followed by £ — cx, bli; and b > by |34, |35); 


3. The recent LHC bounds on the Higgs sector from the measurement of the signal rates 
of YY, TT, WW and ZZ? bd d. 


4. The constraints from B-physics, such as b > sy and B, > pup m 
5. The electroweak precision observables like Mw, sin? 0w and Rẹ " 


In order to reduce the number of free parameters, for the gaugino masses we assume the 
grand unification relation 3M;/5o4 = M2/a2 = M3/a3 and thus we have only one gaugino 
mass parameter (we choose M» in our calculation). In order to explain the muon anomalous 
magnetic moment at 2c level for moderate tan 8 (< 20), for the smuon sector we assume 
the soft-breaking parameters to be 100 GeV bi. For other soft-breaking parameters in the 
squark and slepton sectors, we assume them to be 1 TeV except that for the MSSM we allow 
the third-generation squark mass parameters m; to vary in a wide range. The parameter 
M, in the superpotential of the nMSSM is also fixed to be 1 TeV. Other parameters are 
scanned in the following ranges (among which M, is the mass of the cos Gy, + sin f field, 


i.e., the diagonal element of the doublet in the CP-odd Higgs mass matrix): 


1. For the MSSM: 1 < tan 8 < 20, 100 GeV < u < 600 GeV, 10 GeV < Mə < 200 GeV, 
500 GeV < M, < 3 TeV, 100 GeV < M; < 2 TeV, —3 TeV < A» < 3 TeV. In this 
space we scan two million random points and about three thousands points survived 


the experimental constraints. 


? Note that the constraints from such an invisible Z decay are stringent for a wino-like or higgsino-like 


neutralino, but become quite weak for a bino-like or singlino-like neutralino. 


2. For the NMSSM: 0.1 < A < 0.7, 0.1 < k < 0.5, 1 < tan B < 4, 100 GeV < (u, Ma) < 
1 TeV, 50 GeV < Mə < 150 GeV, 0 < Ay < 1 TeV, —500 GeV < A, < 0. In this 


space we scan ten million random points and about one thousand points survived. 


3. For the nMSSM: 0.1 < A < 0.7, 1 < tan < 10, 50 GeV < (u, Mo, Ma) < 1 TeV, 
0 < Ay <1 TeV, 0 < Mg < 500 GeV, —1 < £p < 1. In this space we scan one billion 


random points and about 2 thousands points survived. 


Note that our scan is not a general scan over the entire parameter space. Since our purpose 
is to figure out if the invisible decay of the Higgs boson is accessible at the LHC, we only 
scanned over a part of parameter space which is potentially able to allow for a light neutralino 
and hence the Higgs can decay into the neutralino pair. 

In Figs. [1] and [2] we display the scatter plots of the parameter space which survive all 
constraints. Fig. []lshows the Higgs decay branching ratios and the LHC (7 TeV) di-photon 
rate, while Fig. 2]shows the observability of the decay h + x?y? through three production 
channels at the LHC (7 TeV). The 2c sensitivity of the ATLAS detector shown in Fig. lis 
obtained from a Monte Carlo simulation carried out in mm For the decay  — xx, the 
signature of the production via vector boson fusion VV — h is two far forward and backward 
tagging jets of moderate pr with considerable missing momentum jy in the central region. 
For the productions pp + hZ and pp — htt, the signatures are obvious: for the former it is 
two isolated high pr leptons from the Z-boson decay and large missing øp from the Higgs 
decay; for the latter it is di-leptons (or lepton plus jets) and large missing pr. 


From these figures we obtain the following findings: 


e In each model the SM-like Higgs can have a mass near 125 GeV, as hinted by the 
recent LHC results. 


e In the MSSM the SM-like Higgs boson dominantly decays to bb (just like in the SM), 
the decay h — x?x? has a very small branching ratio (below about 10%), and the 
di-photon signal rate is close to the SM value. Due to the small branching ratio, the 


Higgs decay h + x9?X? is far below the detectable level. 


e In the NMSSM the decay hı > x9x? can be comparable to hı — bb. In the region with 
a sizable decay ratio of hı + y?y?, the lightest neutralino x? is rather light (below 
h4/2) and the coupling of hı to x? is large (see Eq. [[8). The diphoton signal rate can 
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FIG. 1: The scatter plots of the parameter space which survive all constraints listed in the text. In 
the upper frames, the samples denoted by crosses (sky-blue) show the branching ratio of h + y9x9 
while the samples denoted by circle (magenta) show the branching ratio of h — bb. The solid 
curves (blue) denote the SM prediction for the branching ratio of  — bb. In the lower frames, the 
samples denoted by times (red) show the ratio osysy (pp > h > yy)/osm(pp > h > yy) at the 
LHC (7 TeV). In this figure and below, ’hgjy’ denotes the Higgs boson in the SM, ’h’ denotes the 
lightest nutral Higgs boson in the MSSM, and "^4! denotes the SM-like Higgs boson in the NMSSM 
and nMSSM (the doublet component of hı is over 60%). 


be sizably deviate from the SM prediction, either enhanced or suppressed significantly. 
In a large part of the parameter space, the Higgs decay hı > x?x? is accessible at the 


LHC. 


LHC(7TeV) VV>h 
MSSM 


FIG. 2: Same as Fig.1, but showing the quantity E x Br(hı > y$x9) with esusy (osm) being 
the SUSY (SM) Higgs production rates for the processes VV — h, pp > hZ and pp — htt. The 
solid curves show the 2ø sensitivity of the ATLAS detector at the LHC (7 TeV) with 10 


fbt, 30 fbt and 100 fb~! (the region above each curve is the observable region). 


e In a major part of the parameter space in the nMSSM, the decay hı — y?y? is 
dominant over hı — bb and thus observable at the LHC. The reason is the lightest 
neutralino is singlino-like and is always light, as can be seen from the neutralino mass 


matrix in Eq. (21). Also, from Eq. we see that the coupling gy, oo can be large 
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(near unity). Due to the new sizable decay hı — x?x?, the total width of the SM- 
like Higgs is greatly enlarged and thus its di-photon signal at the LHC is severely 
suppressed. So, if the recently observed di-photon signals at the LHC is verified in the 


near future, this model will be excluded. 


IV. CONCLUSION 


We examined the Higgs decay to neutralino dark matter in low energy SUSY by con- 
sidering three different models: the MSSM, NMSSM and nMSSM. We considered current 
experimental constraints at 2c level (including the muon g — 2 and the dark matter relic 
density) and scanned over the parameter space of each model. Then in the allowed param- 
eter space we calculated the branching ratio of the SM-like Higgs decay to neutralino dark 
matter and examined its observability at the LHC by considering three production channels: 
the weak boson fusion VV — h, the associated production with a Z-boson pp — hZ + X or 
a pair of top quarks pp — htt + X. Our findings are: (i) In the MSSM such a decay is far 
below the detectable level; (ii) In the NMSSM it is accessible in a sizable part of parameter 
space; (iii) in the nMSSM it is detectable in a major part of the parameter space. (iv) When 
this invisible decay is sizable, the visible di-photon decay is suppressed. So, we conclude 
that at the LHC the interplay of detecting such an invisible decay and the visible di-photon 


decay may allow for a discrimination of different SUSY models. 
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